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We study the dynamics of a particle in a horizontaUy and periodically shaken box as a function 
of the box parameters and the coefficient of restitution. For certain parameter values, the particle 
becomes regularly chattered at one of the walls, thereby loosing all its kinetic energy relative to 
that wall. The number of container oscillations between two chattering events depends in a fractal 
manner on the parameters of the system. In contrast to a vertically vibrated particle, for which 
chattering is claimed to be the generic fate, the horizontally shaken particle can become trapped 
on a periodic orbit and follow the period-doubling route to chaos when the coefficient of restitution 
is changed. We also discuss the case of a completely elastic particle, and the influence of friction 
between the particle and the bottom of the container. 
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1. INTRODUCTION AND SUMMARY 



While vertically shaken granular materials have been 
the object of intensive research in the past years, the in- 
vestigation of horizontally shaken granular materials has 
just started Vertically shaken materials show var- 

ious cellular patterns and localized oscillons and 
their horizontal counterpart was recently found to dis- 
play ripple-like patterns [Q. Since these patterns are 
due to the collective behaviour of many interacting par- 
ticles, the one-particle system shows completely different 
phenomena that are, however, equally fascinating. The 
dynamical evolution of a bouncing ball on a vibrating 
platform was studied in As long as the coefficient 

of restitution is smaller than one, a particle that hits the 
platform with sufficiently small relative velocity bounces 
off the platform infinitely often during a finite time and 
looses its memory of earlier dynamics. The authors of |^] 
argue that this "chattering" is the fate of generic trajec- 
tories, which therefore become periodic. They conclude 
that true chaos cannot be observed in this system. 

In this paper, we study the dynamics of a singe particle 
in a horizontally shaken box. While chattering occurs for 
part of the parameter values and initial conditions, we 
find as well other generic scenarios like periodic orbits 
without chatter, the period-doubling route to chaos, and 
strange attractors. The interplay between these different 
modes of behaviour makes this apparently simple system 
astonishingly rich and fascinating. 

The outline of this paper is as follows: In the next 
section we define the system used in our simulations. In 
section ^ we discuss the limiting case of a completely in- 
elastic particle that assumes the velocity of the wall at 
each collision. When the particle hits the wall during the 
half period where the wall accelerates toward it, it sticks 



to the wall until the end of the half period. The num- 
ber of refiections until this locking occurs depends in a 
fractal manner on the parameter of the system. We also 
discuss the influence of friction between the particle and 
the bottom of the container. Next, the opposite limit 
of a completely elastic particle is considered in section 
This system displays all signatures of Hamiltonian 
chaotic systems, including periodic, quasiperiodic, and 
chaotic orbits. The physically most relevant case of a 
partially elastic particle is then studied in section ^ On 
the one hand, increasing the coefficient of restitution from 
zero, the period-doubling route to chaos is observed in 
many cases. Ultimately, the strange attractor becomes so 
large that it includes the chattering region, thus making 
the orbits again periodic. On the other hand, decreas- 
ing the coefficient of restitution from one, the neutrally 
stable fixed points become attractive. All irrational tori 
disappear and all trajectories appear to first become pe- 
riodic. A trajectory starting in the locking region may or 
may not lead back to it, and several periodic orbits coex- 
ist, their basins of attraction being strongly interwoven. 
Section pi concludes the paper. 



2. THE MODEL 

The left and right wall of a horizontally shaken con- 
tainer of length L are described by the equations 

a^icftwaii = Asm{ujt) 
aJrightwaii = L + Asm{ujt) . 

uj is the frequency, and A the amplitude of shaking. We 
denote position and velocity of the particle by x and v 
respectively, and introduce the relative position and ve- 
locity with respect to the container walls 
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Asm{ujt) 



V — Auj cos{ujt) 



(1) 



Between collisions, the particle moves according to a lin- 
ear friction law dv/dt ~ — 7u, i.e., 



du/dt = —"fu + Auj^ sin(a;i) , 



(2) 



with a friction coefficient 7 which will be set equal to zero 
in most parts of this paper. 

The collisions with the wall are partially inelastic, and 
the relative velocity u changes according to 



-rju 



(3) 



at each collision, where rj is the coefficient of restitution. 

It is convenient to measure the particle position in 
units of the amplitude A, time in units of the inverse 
frequency ijJ~^, and velocity in units of Aui. Then the 
system can be described by the dimensionless parame- 
ters a — L/A, 7 = l/^jJ, and 77. We denote the dimen- 
sionless time, length, and velocity again by and u. 
We also introduce the phase of the container oscillation 
(j) = cut mod 2tt. 

The motion of the particle can be written down im- 
mediately. It is most convenient to describe the particle 
dynamics in terms of a map that gives the phase of the 
container oscillation and the particle velocity immedi- 
ately after an impact as function of their values imme- 
diately after the previous impact. Let the particle leave 
the left wall (/ = 0) after an impact at i = <o with a 
relative velocity uq. The particle moves according to 

ko,UQ{t) = 4(costo + uo) - -r^^— (sini -f 7Cost) 
7 7^ -|- i 



1 



Mo 



72 + 1 



(7 sin to — cos to) 



X- exp(-7(t - to)) 
7 



(4) 



the corresponding velocity being u = dl/dt. We will 
mainly discuss 7 = 0, where this equation reduces to 

lto,uo{t) = sin to - sint -I- {uq + costo)(t - to) . (5) 

The next collision occurs at the smallest solution tc > 
to of 



impact at same wall (6a) 
impact at other wall (6b) 



at which the velocity is Uc — utg^uoitc)- If the particle 
impacts again the same wall, the new velocity is now 
Ml = —tjUc- On the other hand, if the particle collides 
with the right wall, it is convenient to use the symmetry 
between left and right wall to map the right wall to the 
left wall via t —^ t—ir and I —^ a— I, so that now ui — rjUc- 
The velocity is thus always measured with respect to the 
wall at which the last impact has occurred. 



Given an initial phase ( 
arrive at the map 

d>-[ = tr mod 27r . 



to and velocity uq, we thus 



ui = -rjUc 
impact at same wall (7a) 
01 = {tc — tt) mod 27r , ui = rjUc 

impact at other wall (7b) 

Before we will study this mapping in detail below, we 
shall briefly discuss the case of no dissipation (77 = 1 
and 7 = 0). Here, the dynamics_can be derived from a 
time-dependent Hamiltonian 

H ^ j- + VsQ{q)+AmLu\sm{Lut), (8) 

where Vsgiq) is a square-well potential describing the 
walls at g = and L, and m and p are mass and momen- 
tum of the particle. The last term in the Hamiltonian (||) 
causes a periodic oscillation of the particle with respect 
to the walls. (Since we describe the system in terms of 
relative coordinates and velocities between the particle 
and the wall, a system with periodically oscillating walls 
and a system with a periodic particle oscillation superim- 
posed to the ballistic motion are equivalent.) Measuring 
p in units ofmLuj, H and V in units of toL^w^, q in units 



of L, and t in units of w ^, 
formulation. 



we arrive at a dimensionless 



H 



p^ 
2 



'->SQ 



(9) 



■ sint 



(9) 



where vsqiq) is an appropriately scaled square well po- 
tential with walls at and 1. Upon introducing action- 
angle variables (J, 9) for the "unperturbed" Hamiltonian 
^^0 = f^ + vsgiq) the Hamiltonian (||) transforms 
to 



H = 



(1 2 °° 1 ^ 

2^'"*+^ E _sin(n0-t) 



(10) 

The transformation between (p, q) and (J, 9) is given by 

p = 7rJsign0 and q = \9\/Tr (11) 

for — TT < 9 < TT, and periodically continued in 9. We 
have included the term (2a) ~^ sint in the Hamiltonian 
( pO| ) for completeness sake, but we shall drop it below as 
it does not influence the dynamics. 

3. THE COMPLETELY INELASTIC PARTICLE 

3.1. Modelling by a one— dimensional map 

We now discuss the case 7 = and r/ = in detail. 
The particle moves freely between the walls, and after an 
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impact the relative velocity uq with respect to the wall is 
zero. This means that the two-dimensional mapping de- 
fined above is in fact reduced to a one-dimensional map- 
ping. The subsequent fate of the particle depends only 
on the phase 0o at the moment of impact, as indicated 
in Fig. g. 




absorbing region 

< > 



FIG. 1. The fate of a particle that is inelastically reflected 
at the wall. 

For e [0, 7r/2[ mod 2tt, the particle is reflected from 
the wall with (absolute) velocity cos((/)o) and is headed 
for the other wall. For <j)o e [k,2tt[ mod 2tt, the particle 
sticks to the wall until (po — 2tt, and then leaves the wall 
with (absolute) velocity 1. In the intermediate region 
00 €:]7r/2,7r[ mod 2tt, the sign of the particle velocity is 
not reversed during the collision, and the particle hits the 
same wall again at a later phase which is determined 
by the equations (||), (|6a|), and ( |7a]) , leading to 

sin 00 — sin 01 -I- (01 — 00 ) cos 00 = . (12) 

For 00 > 0^ = 1.79.., the phase 0i is in the locking inter- 
val [tt, 27r[. As we are interested in the dynamics between 
impacts at opposite walls, we can extend the locking in- 
terval to include 0^, i.e., to the interval [0c,27r]. For 
00 < 0°, the phase 0i is in the interval [0,7r/2[ mod 2tt 
and the particle is immediately reflected towards the 
other wall, where it arrives at a phase 02 now given by 
(6b) and (7b), i.e., 

sin 01 + sin 02 4- (02 — tt — 0i) cos0i = a . (13) 

From now on, let 0^ denote the phase of a wall at the 
moment where the particle is reflected towards the other 
wall for the ith time. We do not count reflections that 
lead to a subsequent reflection at the same wall. Then, 
the map that gives 0i+i as function of 0^ is a map of the 
interval [0, 7r/2[ into itself, dependent only on the dimen- 
sionless parameter a — L/A. The map 0^+1 — $Q,(0i) 
is given implicitly by the equations (|lj) and (|l^), and is 
shown in Fig. |^ for the value a = 5. Branches with pos- 
itive slope correspond to a particle arriving at the other 
wall at a phase in the interval [0,7r/2], from where it is 
immediately reflected back, while branches with negative 



slope indicate a twofold reflection at the other wall, the 
first reflection being in the phase interval ]tt/2, 0^^]. The 
dashed line is the diagonal 0i+i — (f>i, and the dotted line 
indicates the fate of a particle starting at phase zero. For 
tt/2 > 0,; > tt — 0^^, the velocity of the particle is so small 
that it cannot hit the other wall in the locking region. A 
small interval A0i maps onto an interval A0i+i propor- 
tional to the flight time, which is of the order a/cos0i. 
For this reason, the density of branches in the map in- 
creases with increasing a and 0, and diverges when 0^ 
approaches the value 7r/2, where the initial velocity de- 
creases to zero. 




0.0 0.5 1.0 1.5 



FIG. 2. The map for a = 5. 

3.2. Periodic trajectories 

The phase $0,(0) at which a particle starting at phase 
zero will be reflected from the other wall depends on a. 
It is zero for a S [{2n — 1)tt + (p'^ + sin0^, {2n + 1)tt], for 
any nonnegative integer n. In this particle that 

starts in the locking region will hit the locking region at 
the other wall, and the periodicity of such a trajectory is 
one. For all values of a, the map intersects the perpendic- 
ular axis with a slope zero, since the initial velocity does 
not change with 0^ for 0^ = 0. For values a only slightly 
above (2n -I- l)7r, the map therefore has a stable fixed 
point (f)* close to zero. This fixed point vanishes with in- 
creasing a via a saddle-node bifurcation, and for values 
of a slightly beyond the bifurcation, a particle trajectory 
can be trapped for a long time in the neighborhood of the 
former fixed point, before it escapes and hits ultimately 
the locking region where 0i+i — 0. Due to these lock- 
ing regions, there exist no truly chaotic trajectories, but 
each trajectory has a finite periodicity. When a increases 
further, the number of reflections in the neighborhood of 
the former fixed point (or, more precisely, on the first 
branch that has positive slope) , decreases in steps of size 
one. Close to such a decrease, a trajectory that leaves the 
first branch goes through the upper right-hand corner of 
the map and can therefore have an arbitrarily large pe- 
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riodicity. Further away from it, a trajectory can hit the 
locking region immediately after leaving the first branch. 
In general, very complex trajectories can occur. In par- 
ticular when ^a{0) is close to 7r/2, the trajectory will 
spend some time in the upper right-hand corner of the 
map, and a slight change in a may lead to a large change 
in the periodicity of the trajectory. 



There are also many possibilities for obtaining low- 
periodicity trajectories. In Fig. ||, e.g., the periodicity 
is two. A period-doubling bifurcation sequence occurs 
when a is decreased from (2n — 1)tt + 0'^ + sin^'^. In 
Fig. ^, this scenario is sketched. Since it involves only 
the first two branches of the map, the other branches 
are not shown. As long as the first branch is small (e.g., 
for a = 5.8), the particle hits the locking region after 
the second reflection. With decreasing a, the endpoint 
of the trajectory moves to the right, and finally hits the 
foot of the second branch, leading to a period-doubling 
bifurcation. For a = 5.743, one observes therefore a cycle 
of period 4. When a is decreased further, the end point 
of the 4-cycle moves to the left and hits the foot of the 
first branch, leading to another period doubling. 

1.5 I • , • , 
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FIG. 3. A series of bifurcations when a is decreased. 

The last part of Fig. shows an 8-cycle for a = 5.7406. 
This period-doubling scenario continues as the end point 
of the 8-cycle moves to the left, etc., and the period be- 
comes infinite when the trajectory hits the unstable fixed 
point on the first branch. For even smaller values of a, 
other periodicities are observed that are not powers of 2. 

Fig. ^ shows the period of a trajectory originally start- 
ing in the absorbing region as function of a. 
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FIG. 4. Period of trajectories originally starting in the ab- 
sorbing region as function of a. 

All the above-mentioned features can be seen: The 
large plateaus of period 1 that extend beyond (2n-|- 1) x tt 
because of the existence of stable periodic orbits outside 
the locking region; the stepwise decrease of the num- 
ber of subsequent reflections at the first branch as a 
increases, and a high periodicity in between these sub- 
sequent plateaus; the trajectory of period 2 shown in 
Fig. 1^; the period-doubling sequence at the left end of the 
period-1 plateaus (due to the finite resolution in a only 
the 2- and 4-cycles can be seen. The 3-cycle lies beyond 
the period-doubling sequence). As the branches of the 
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map become steeper and more numerous with increasing 
a, plateaus will become shorter, and the structures will 
become more complex for large a. 



3.3. Comparison to a vertically shaken particle 



In |g|7|, the periodicity of trajectories of vertically 
shaken completely inelastic particles is studied. There, 
gravitation brings the particle back to the vibrating plat- 
form, and no second wall is needed. Therefore, fast par- 
ticles have longer flight times between reflections than 
slow particles, while the opposite holds for horizontally 
shaken particles. For large oscillation frequencies, the 
authors of obtain an approximate map of the form 
~ -f const. X cos An approximate map for 
our system is obtained from Eq. (|l^) for large a and 
reads (pi+i — (j>i + tt + a/ cos (j>i. The map in ||^ has 
no stable fixed points and no region where the trajec- 
tory can be temporarily trapped. Nevertheless, there are 
infinite hierarchies of periodic trajectories between low- 
period plateaus, leading to a structure on all scales of the 
period-vs-frequency plot, as for the horizontally shaken 
particle. 



3.4. The influence of friction 

In the presence of friction between the particle and the 
container bottom, there exist stable low-velocity particle 
trajectories that never touch the walls, but move period- 
ically back and forth in a small region somewhere in the 
central part of the container. Particles that are refiected 
at a phase near 7r/2 have a very small velocity and can 
become trapped in such a region. Apart from the locking 
sections in the map Fig. |^, there is consequently an ab- 
sorbing section for (f)i above some threshold that depends 
on a and the strength of friction. From equation (^) with 
Mo = we see that the maximum possible value of I is 



00,0 



— COS 00 

7 



1 



(14) 



The first term is the initial velocity divided by the friction 
coefficient, and is the distance that a particle with initial 
velocity cos((/)o) can travel on a stationary surface. The 
second term is the amplitude of the periodic oscillation 
that a particle performs when it is periodically driven. 
If Z^j^Q < a, the particle never reaches the other wall. 
If a > 1, this is always the case for (/)o close enough to 
7r/2, as the initial (absolute) velocity of the particle is 
arbitrarily small. If a is large enough, i.e.. 



a > 



1 



1 



7 yiT^ ' 



(15) 



then even a particle that was reflected at phase zero, i.e., 
with maximum velocity, does not reach the other wall. 



Fig. H shows the period of trajectories as function of a 
for 7 = 0.025. 
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FIG. 5. Period of trajectories originally starting in the ab- 
sorbing region as function of a, with a friction coefficient 
7 = 0.025. 

Trajectories that get stuck in the absorbing region are 
assigned a period zero. To keep the computer program 
simple, friction was only active for x-coordinates between 
1 and a — 1, and only leading terms in 7 were considered. 
Comparing Fig. ^ with Fig. ^ one finds that with increas- 
ing a the graph becomes more and more compressed hor- 
izontally, and that an increasing fraction of trajectories 
get stuck. The region a > 41.0 where all trajectories get 
stuck, is not shown on the graph. 



4. THE COMPLETELY ELASTIC PARTICLE 

We now study the non-dissipative case rj — 1 and 
7 = 0. Here, the particle moves without friction and 
is reflected elastically at the container walls, where its 
relative velocity with respect to the wall changes sign at 
each collision. The velocity is now an independent vari- 
able, adding a new dimension to the phase space. As the 
system is non-dissipative, it can be described by a time- 
dependent Hamiltonian as described in section An 
equivalent problem is studied in |p|-p^, namely the mo- 
tion of a charged particle in a one-dimensional infinite 
square well potential driven by an oscillating external 
field. In those papers, the particle performs periodic os- 
cillations (superimposed to the free motion) . The param- 
eter a^^ is proportional to the electric field. In terms of 
coordinates and velocities relative to the container walls, 
the dynamics is the same as for our system. 

Since the system is hamiltonian, phase space volume 
is conserved in time, and a stroboscopic Poincare section 
that gives l{t) and u{t) in time intervals T = 27r is area 
preserving. However, we find it physically more relevant 
to plot the phase of the oscillation and the velocity of 
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the particle relative to the wall immediately after each 
reflection, i.e. (j) = tc mod 2tt and u{tc) for l(tc) = 0. As 
the time intervals are not constant, this representation 
is not area preserving, but describes of course the same 
physics as a true Poincare section. Fig. ^ shows our data 
points for a — 1.65, a — 16.5, and a = 165, and for sev- 
eral different trajectories. The value of a used in |^-[Tc|] 
is of the order of several hundred, since the authors were 
mainly interested in the limit of small electrical field. 

In order to understand some important features of 
Fig. 1^, we turn to the action-angle variable representation 
® , i.e. 



H = 



J2 -^sin(n^-i) 



If one now attempts a canonical transformation such 
that the new action and the new Hamiltonian arc con- 
stants of motion to first order in a"^, one finds that this 
perturbation theory breaks down near the resonances, 
where n^Jn = 1, i.e., u = a/im, for any positive n. For 
large a, the resonances to lowest n (i.e., to highest en- 
ergy), are widely separated, and in the neighborhood of 
such a resonance, the system is well approximated by a 
single resonance Hamiltonian 



Hn 



-K^J^ 2 sin(n0 - t) 



(16) 



This Hamiltonian is integrable and becomes identical 
to the Hamiltonian for a pendulum under a canonical 
transformation that equals 6 — t with the pendulum an- 
gle variable. For a — 165, Fig. shows clearly these 
pendulum-like trajectories. In Pi, it was shown that 
the breakdown of the irrational tori between two reso- 
nance regions and the onset of large-scale chaos can be 
approximately determined from a two-resonance Hamil- 
tonian that includes the two resonances adjacent to the 
considered torus. This chaotic region can be seen in the 
lower part of all three Figs. |^. The stable periodic orbits 
that can be seen for the larger two values of a represent 
trajectories that hit the walls always at the moment of 
farthest extension, as illustrated in Fig. ^ and they are 
at the center of resonance zones and correspond (in the 
single resonance approximation) to a pendulum at stable 
rest. The velocity of such a trajectory is given by 



(17) 
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FIG. 6. Phase and velocity relative to the wall at the mo- 
ment of reflection for a — 165, 16.5, and 1.65, for several 
different trajectories. To each trajectory, several 100 points 
are plotted. 



The stability of these trajectories is found from linear 
stability analysis. The eigenvalues of the stability matrix 
are 



ui V ul ul 



(18) 
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FIG. 7. The periodic orbit of lowest period. 

For > 7r/2 the eigenvalues are complex conjugate 
with |A| = 1, and the trajectory is neutrally stable. For 
< 7r/2, one eigenvalue is in modulus larger than one, 
and the trajectory is unstable. To each of these stable 
periodic orbits there exists an unstable counterpart that 
hits the walls always at phase 7r/2 and that has the veloc- 
ity wj^ = (a — 2)/(2ri— l)7r. These trajectories correspond 
(in the single resonance approximation) to a pendulum at 
its unstable fixed point. For nonintegrable systems, un- 
stable periodic orbits are always part of chaotic regions. 
It is, however, remarkable, that the boundary of the large 
chaotic regions is close to such an unstable fixed point, 
and that its shape is close to the shape of a pendulum sep- 
aratrix u ~ Ci -I- C2VI ~ sin <j). This means that the sep- 
aratrix belonging to the highest-energy resonance within 
the chaotic region determines the approximate boundary 
of the chaotic region. 

For very large velocities u ^ a, the phase and veloc- 
ity relative to the wall change only by a small fraction 
between subsequent reflections, and solving a differential 
equation to leading order in l/u, one finds on approxi- 
mate equation for the tori 



u ~ \/ C — a sin . 

(More generally, for a container oscillating according to 
an equation ajwaii = /(</>) i one finds u ~ yj C — af"{4}).) 
This approximation can also be found from the Hamilto- 
nian (^0|) . A more detailed calculation including the first 
resonance gives 



it~^^^^ ^l + C"±yC"- ^sin</>j -asin(/) (19) 

The shape of stable irrational tori for smaller values of 
the velocity is obtained in pO| ] through approximation by 
periodic trajectories. 



5. THE PARTIALLY INELASTIC PARTICLE 



The dynamics of the partially elastic particle (0 < 77 < 
1) show a very intricate structure. Depending on the pa- 
rameters and the initial state, the trajectory of a particle 
can be on a low-periodicity periodic orbit, a strange at- 
tractor, or an orbit that looses periodically all its energy 
due to "chattering" , i.e. an infinite number of reflections 
at the same wall during a finite time interval. In the fol- 
lowing subsections we discuss all these phenomena and 
the interplay between them as the parameter values are 
varied. 



5.1. Chattering 

When the particle hits the wall with sufficiently small 
relative velocity, it is reflected infinitely often from that 
wall during a finite time interval and looses all its energy, 
as described in [Q for the vertically shaken particle. It 
subsequently sticks to the wall until the phase is 4> — 0, 
and then leaves the wall with relative velocity u = 0. 
In order to understand this phenomenon, let us consider 
first a particle colliding with a wall that is accelerated at 
a constant rate a. This situation is equivalent to a ball 
bouncing in a gravitational field, where it experiences a 
constant acceleration toward the ground. If the initial 
relative velocity of the particle is mq, the relative velocity 
becomes zero after the time 

T = 2rjua/ a + 2ri'^uo/a + .... = 2uo7y/a(l — 77) 

after the first collision. In our system, the wall accelera- 
tion is of the order Au'^, and chattering can occur if the 
time T is no larger than of the order of the oscillation pe- 
riod u>~^, leading to the condition uq ^ Auj{l — rj)/rj for 
chattering. For 77 close to 1, uq must be very small, which 
is only possible if the particle hits the wall at a phase (/)o 
within a distance A(/)o oc ^1 — rj of tt, or, equivalently, 
if the particle's trajectory hits the (absolute) position 
X = (or X — L) within a time interval of the order 
Ai ~ uoA(/)o cx (1 — 77)^^^. The phase space volume 
1*0 Af for which chattering occurs shrinks consequently 
as (1 — 77)^/^, when 77 approaches 1. 

In 0, it is argued that chattering is the generic fate 
of a particle bouncing on a vibrating platform, with the 
number of reflections between two chattering events di- 
verging as rj approaches 1. In our system, however, a 
particle that leaves the wall at phase ^ = and with 
relative velocity u = 0, returns to the chattering region 
only for certain combinations of a and 77, and is oth- 
erwise trapped on a periodic or chaotic orbit. Even in 
those cases where a chattered particle is chattered again, 
there may exist other trajectories that never enter the 
chattering region. 
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5.2. The almost elastic case 



Due to dissipation, a given phase space volume shrinks 
with time. In particular, regions with large initial veloc- 
ity are depleted, since the energy of a particle decreases 
exponentially fast until it reaches the regime where the 
time between impacts is of order unity, i.e. where its ve- 
locity is of order a. The change of u with the impact 
number n is approximately given by 

du/dn = — u(l — 77) . 

Using the relation dcf) ~ {a/u)dn, we find 



(20) 



for the decrease of u with 4>, with a constant C that 
is determined by the initial conditions. A small initial 
velocity interval consequently increases rapidly with in- 
creasing (/i, leading to a strong shear of an initial domain 
of phase space, and interweaving different domains of at- 
traction, as we shall see below. In fact, every domain of 
attraction contains points with arbitrarily large velocity. 

Another consequence of the shrinking phase space vol- 
ume is that all stable periodic orbits of a Hamiltonian 
system become attractive when small dissipation is added 
[ pl| . This is because the eigenvalues of the stability 
matrix remain complex conjugate for sufficiently small 
(1 — 77), while their product is A1A2 = 77^, leading to 
I Ail — r] < 1. Thus, for rj infinitesimally smaller than 
1, there is a plethora of attracting periodic orbits with 
competing domains of attraction. When ry decreases fur- 
ther, most of these periodic orbits merge with unstable 
orbits and disappear, leaving only few attractive fixed 
points. The fixed points in the center of the resonances 
survive over a fairly large range of rj values, as shown by 
the subsequent calculation. Writing 



p={2n- l)7r 



1 



1-7J 



(21) 



the fixed point in the center of the nth resonance changes 
with decreasing rj according to 



sm = 7i and u = cos ( 



(22) 

This fixed point can vanish by merging either with an 
unstable fixed point, or it can become unstable. The 
merging of stable and unstable fixed points occurs when 
the square root in (^2|) becomes imaginary, leading to the 
condition 



ri > 



y/a^ - 4 - {2n - l)n 
\/oP- -4-1- (2n - l)7r 



(23) 



for the existence of the fixed point. 
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FIG. 8. Pliase and velocity relative to tlie wall at the mo- 
ment of reflection for a — 165 and = 0.9999, for three 
different trajectories. To each trajectory, 50000 points are 
plotted. 

On the other hand, a stability analysis gives a stable 
fixed point for 



?7 



1-77 (271 - l)n 



< tan (h < 



1 + 77 (2ri - l)7r 



(24) 
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The violation of the left-hand side of condition (24) cor- 
responds to a period doubling bifurcation (one eigenvalue 
becomes —1), the violation of the right-hand side corre- 
sponds to a saddle-node bifurcation (one eigenvalue be- 
comes 

The chaotic regions surrounding the unstable fixed 
points for f] — I also continue to influence the dynamics 
for 77 < 1. Typically chaotic regions become tran- 
sient when a small dissipation is added, and their main 
effect is the mixing of basins of attraction. In figure ^ we 
illustrate this for 77 ~ 0.9999, where we show three dif- 
ferent orbits of particles starting at high initial velocities 
in the case a = 165. One clearly sees the intertwining of 
different domains of attraction above the first resonance, 
making the fate of each particle appear random. After a 
long transient time, the first trajectory becomes trapped 
in the first resonance, the third trajectory in the second 
resonance, and the second trajectory is not yet trapped 
after 50000 refiections. It might either become attracted 
later by a higher-period orbit, or it might be chattered 
regularly. In the latter case, chatter should occur at inter- 
vals of the order of (1 — fy)"^^^ ~ 3 x 10^^ refiections (see 
previous subsection). A region that is chaotic for ~ 1 
is not immediately transformed into a strange attractor 
when r] is decreased However, strange attractors 

occur for smaller values of 77, as we shall see below. 



gion, or to a stable periodic orbit that formerly coexisted 
with the attractor. Alternatively, as in the Feigenbaum 
scenario, an increase in 77 can induce a transition from 
a strange attractor to an intermittent periodic regime of 
some odd periodicity (e.g. 3), that undergoes in turn 
a period-doubling scenario. Even in cases where the 
attractor is destroyed due to chattering for some value 
of 7/, such intermittent periodic regimes and subsequent 
period-doubling cascades can still be observed when 77 is 
further increased. 
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5.3. Large dissipation 

As we have seen in Section |[ every trajectory is peri- 
odic for 77 = and may or may not go through the chat- 
tering region (f> e [tt, 2tt[ mod 2tt. Both types of orbits 
can coexist for small 7/ if the trajectory that starts in the 
chattering region is not trapped by the other periodic or- 
bit. Since the basins of attraction of the chattering orbit 
and the normal periodic orbit are strongly interwoven, a 
small change in 77 can induce or remove such a trapping. 
When 77 is increased, normal periodic orbits typically go 
through a period-doubling scenario and then become a 
strange attractor. This attractor may be a global attrac- 
tor, or it may coexist with other attractors, e.g., periodic 
orbits. Due to the large dissipation, the attractors are 
rather fiat and seem almost one-dimensional in Poincare 
like plots, as also known from other systems with large 
dissipation, like the Lorenz attractor [g2j . Fig. ^ shows 
a strange attractor for a = 10 and 77 — 0.142. The lower 
arches of the attractor correspond to the second, third, 
and fourth refiection at the same wall and are not present 
for smaller rj. When rj increases further, the number of 
arches diverges, and the attractor is replaced by a chat- 
tering orbit. 

There are other mechanisms by which strange attrac- 
tors may vanish with increasing rj. As the attractor in- 
creases with rj, its edge points may approach an unstable 
periodic orbit that ultimately diverts the trajectory from 
the attractor and leads it either into the chattering re- 



FIG. 9. The strange attractor for a = 10 and rj — 0.142. 

Figure ^ shows the various attractors for a = 2 as 
?7 is increased from to 1. Increasing 77 from zero, we 
see a stable fixed point emerging from the chattering do- 
main which dominates the dynamics at 77 = 0. This fixed 
point undergoes a complete period doubling scenario as 77 
increases further. At the end of the period doubling sce- 
nario we have a chaotic attractor which continues to grow 
until it disappears when it overlaps with the chattering 
domain. The latter part of the period doubling scenario 
exists simultaneously with the stable three-cycle that is 
also shown in the figure. This three-cycle moves to the 
left as 77 increases, and it survives to 77 = 1. Close to 
77 = 1, more periodic orbits exist that are not shown in 
the figure. Moreover, there are other small attractors 
which appear for intermediate values of rj (some of the 
"dust" in the figure). 

We have chosen a = 2 for Fig. ^ as here the pe- 
riod doubling is clearly visible. For other values of a, 
more complicated scenarios can occur, as already indi- 
cated in the text above. There may be several period 
doubling cascades, which can occur for increasing as well 
as decreasing values of rj. Period-doubling cascades for 
increasing 77 are usually initiated by an eigenvalue of the 
stability matrix of a periodic orbit becoming equal to 
— 1, while cascades for decreasing rj occur when a peri- 
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odic trajectory starts to experience several reflections at 
the same wall. Moreover, there can be intervals of rj for 
which all trajectories experience chattering, i.e., where 
the only attractor is the periodic orbit with chattering. 
Period-doubling cascades are also observed when a is 
varied for fixed rj. 



1.5 




1-6 3,6 5.6 7.6 

phase at reflectioii 

FIG. 10. The set of attractors for a = 2 as ?7 changes from 
to 1. 



5.4. Ghost attractors and domains of attraction 

Since the size of strange attractors increases with in- 
creasing 7y until they become absorbed by the chattering 
orbit, strange attractors become rare for higher values 
of 77. Nevertheless, the "ghosts" of these attractors still 
continue to influence the dynamics by trapping trajecto- 
ries for a transient time, which can last several hundred 
or thousand reflections. Consequently, the domains of 
attraction for normal periodic and chattering orbits be- 
come strongly interwoven, as shown in Fig. |ll| for a = 10 
and ry = 0.61575. There is a fixed point sitting in the 
center of the curl, and there is one chattering orbit. The 
large black and white areas around the curl are the do- 
main of attraction of the fixed point, and the two colours 
distinguish between trajectories that need an even or odd 
number of refiections before reaching some small region 
around the fixed point. The area with the irregularly 
distributed black dots is the region where the domains of 
attraction of the fixed point and the chattering orbit are 
interwoven, the black dots being points that utimately 
become attracted to the fixed point. For two nearby 
starting points in this region, the trajectories follow a 
"ghost attractor", by which they are trapped for quite 
some time. The release of the trajectories to either the 
attracting fixed point or the chattering domain is seem- 
ingly random due to the influence of the attractor, as can 
be seen from the bottom flgure. The white region at the 
bottom of the top flgure belongs entirely to the basin of 
attraction of the chattering orbit. 




phase at mfleotio?^ 
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FIG. 11. The two domains of attraction for a — 10 and 
ri = 0.61578. The attracting fixed point is at <^ ~ 5.7 and 
u ~ 2.7. Chattering occurs immediately for orbits entering 
the white domain at the bottom of the upper figure, and tra- 
jectories starting in the large black and white areas around the 
fixed point become absorbed in a small neighborhood around 
the fixed point after an even resp. odd number of steps. In 
the region with the irregularly placed black dots, trajectories 
are strongly mixed, and it cannot be predicted whether they 
will be chattered or absorbed at the fixed point. The lower 
figure illustrates how how two nearby trajectories follow a 
"ghost attractor" , before they eventually either chatter or get 
attracted to the fixed point. 



5.5. The influence of friction 

For the case of a completely inelastic particle, we have 
studied in Section |^ the influence of friction between 
the particle and the container bottom. Particles that 
arc so slow that they cannot reach the other wall be- 
come trapped in the region between the walls, where they 
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perform a low-amplitude periodic oscillation. For suffi- 
ciently large container sizes, this is the fate of all particle 
trajectories. A similar behaviour can be expected when 
the coefficient of restitution does not vanish. All orbits 
that extend to a region of sufficiently low velocity will 
become trapped in the region between the walls. This 
will in particular affect long orbits like chattering orbits 
and strange attractors. For the completely elastic case, 
even a small amount of friction destroys the hamiltonian 
nature of the system, and we expect that the system 
behaves similar to the situation where the coefficient of 
restitution is slightly smaller than 1. Thus, we expect 
all stable periodic orbits to become attractive, and all 
chaotic regions to become transient for an infinitesimally 
small amount of friction. 



moving wall is studied for approximately ten particles. 
The authors find that most of the particles form a clus- 
ter almost at rest, while a few remaining particles travel 
between the boundaries at a much higher speed. We have 
seen a slightly different phenomenon in the periodically 
shaken box, namely the formation of two clusters travel- 
ling between the boundaries and the center of the system, 
similar to Newton's cradle. 
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6. CONCLUSION AND OUTLOOK 



In this paper, we have shown that a particle in a hori- 
zontally shaken box shows a much richer behaviour than 
a particle on a vibrating platform. While chattering, 
i.e., the loss of all kinetic energy during a finite time, 
may occur for any value of the coefficient of restitution 
smaller than one, other scenarios like period-doubling 
and strange attractors are observed as well. We have also 
discussed the interplay and transitions between these sce- 
narios, and the influence of friction between the particle 
and the container bottom. We know only of one other 
system where the interplay between chattering, chaos, 
and periodic orbits has been studied before, namely the 
forced oscillator impacting on a wall . 

The chattering phenomenon discussed in this paper is 
a good approximation to reality only when the collision 
time between the particle and the wall is much shorter 
than the oscillation period. If this condition is not satis- 
fied, the deformation of the particle during the collision 
has to be included in the model, thus adding another 
degree of freedom. 

When N > 1 particles are placed in the container, new 
phenomena will arise. Three or more particles with a suf- 
ficiently small coefficient of restitution rj < 1 — 1/N are 
known to undergo an inelastic collapse, where they loose 
all their relative kinetic energy due to infinitely many col- 
lisions during a finite time , a phenomenon similar 
to chatter. But even for parameter values that do not 
allow for an inelastic collapse, clustering phenomena oc- 
cur. In ijl^, a system with one elastic and one thermally 
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